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Abstract 

With the recent observation of nonzero #13, five neutrino oscillation parameters are 
now known. By imposing four zeros in the Yukawa coupling matrix of the type I seesaw 
model, the number of parameters in the neutrino mass matrix is reduced to seven, and 
we are able to make predictions for the lightest neutrino mass, Dirac CP phase, and 
neutrinoless double beta decay. Four texture zeros in the Yukawa coupling matrix is 
equivalent to either a single texture zero or a single cofactor zero for an off-diagonal 
element of the light neutrino mass matrix. We find strong similarities between single 
texture zero models with one mass ordering and single cofactor zero models with the 
opposite mass ordering. In the context of a specific class of single-flavor leptogenesis 
models, we find additional constraints on the parameter space. 
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1 Introduction 



Low-energy neutrino phenomenology is described by nine parameters in the Majorana mass 
matrix of light neutrinos, which can be written as [fj M = V*dia,g(mi, m^)V\ with 
V = f/diag(l,e^ 2 / 2 ,e^ 3 / 2 ) and 



U 



C13C12 

— S12C23 — Ci2S23Si3e J 



■512-523 ~ Cl2C23-5l3 e 



/5 



C13S12 Si3 e 

iS 

C12C23 — Sl2S23 s 13 e S23C13 



"C12S23 ~~ 5120235136 C23C13 



(1) 



After the measurement of #13 by the Daya Bay [2], RENO [3], and Double Chooz |3] exper- 
iments, five of them are known; the result of a recent global three-neutrino fit [5] is shown 
in Table 1. The tiny masses of light neutrinos can be elegantly explained by the seesaw 
mechanism [6], but, unfortunately, with the introduction of additional free parameters that 
cannot be measured in the forseeable future. The most popular seesaw model is the type I 
seesaw, in which the mass matrix of the light neutrinos can be written as 



M = X T M R - 1 Xv 2 = Y T Y . 



(2) 



where, for iV heavy right-handed neutrinos, Mr = diag(M 1; M 2 , • • ■ , Mjy), A is a iV x 3 
Yukawa coupling matrix, and 



Y = vM R - 1/2 X 



yie yi(t yir 

V2e V2n V2r 



(3) 



VNe VNp. VNt 

Note that on permuting the rows of the Y matrix (which is equivalent to reordering the right- 
handed neutrinos) or applying a rotation to the rows of the Y matrix (which is equivalent to 
a rotation in the space of the right-handed neutrinos), the mass matrix of the light neutrinos 
remains the same. 

The standard seesaw model has iV = 3 and permits any set of low-energy neutrino pa- 
rameters in M u . A way to extract predictions from the seesaw model is to impose constraints 
on its parameters. The most economical seesaw model includes two right-handed neutrinos 
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Table 1: Best-fit values and 2a ranges of the oscillation parameters [5], with 5m 2 = m\ — m\ 
and Am 2 = m 2 — (m\ + m|)/2. 



Ordering 


M°) 


M°) 


M°) 


5m 2 (10" 5 eV 2 ) 


|Am 2 |(10- 3 eV 2 ) 


Normal 


33.6l 2 2 ;J 


8.9^| 


38.41^ 


7 ^4+0.46 
' - d ^-0.39 


2 4S+ - 12 


Inverted 


33.6l|J 


9.0i?;« 


38.810 ©47.5 - 53.2 


7 e;4+0.46 
'•°^-0.39 


2 42+ - 11 



and two zeros in the Yukawa coupling matrix, or, equivalently, two zeros in Y [7j. The 
measurement of #13 excludes the normal mass ordering (NO) and stringently constrains the 
allowed parameter space for the inverted mass ordering (10); in particular, the phase 8 must 
be such that Dirac CP violation is close to maximal [8J. Also, for two right-handed neutrinos, 
one of the light neutrinos must have vanishing mass, so the allowed values of the sum of all 
light neutrino masses (which affects structure formation in our universe) and \M ee \ (which 
determines the rate for neutrinoless double-beta decay (Ou/3/3), a signal of lepton number 
violation) are limited. 

In this paper, we extend this most economical model to include a third right-handed 
neutrino. We use texture zeros in Y as extra constraints, which may arise in, e.g., extra 
dimensional models [8]. If there are 5 or more texture zeros in Y, the most economical model 
with two right-handed neutrinos or a block diagonal matrix is obtained; the latter is excluded 
by the current experimental data. So the simplest case for three right-handed neutrinos has 
four texture zeros in Y, which is equivalent to five nonzero elements; for previous work 
see Refs. P, [10]. Here we derive analytic formulas that relate the free parameters to the 
dependent ones and determine the constraints on these models including the recent data on 

#13- 

With five nonzero elements in Y only two phases are physical, so that there are seven free 
parameters in the four texture zero model. Hence, we can use the five observed oscillation 
parameters from the global fit to determine the allowed regions for the Dirac CP phase 5 and 
mi (7B3) for the normal (inverted) ordering. Then we can obtain the values of the Majorana 
phases 02 and 03, completing our knowledge of all the elements in the light neutrino mass 
matrix. Furthermore, from Eq. (J2J) we can find Y and then study neutrinoless double beta 
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decay and leptogenesis. 

In Sec. 2, we list the general properties of the four texture zero model and give a 
brief review of neutrinoless double beta decay and leptogenesis. In Sec. 3, we use current 
experimental data to study the allowed parameter regions for each texture. We conclude in 
Sec. 4. 

2 Classes with four texture zeros 

There are 3 basic ways to have 5 nonzero elements: (2, 2, 1), (3, 1, 1) or (3, 2, 0), where the 
numbers indicate how many nonzero elements are in each row of the Y matrix. The (3, 2, 
0) case is equivalent to two right-handed neutrinos with one texture zero element and one 
vanishing mass, and is equivalent to the most economical N = 2 model after a rotation in the 
right-handed neutrino space. So only the (2, 2, 1) and (3, 1, 1) cases need to be considered. 
There are 9 independent textures in the (2, 2, 1) case and 3 independent textures in the (3, 
1,1) case. We divide them into four classes: 
Class 1 
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X 
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X 







X 


X 





1A : 


X 
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IB : 


X 
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1C : 


X 





X 




X 













X 













X 



Class 2 
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X 







X 


X 
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X 





2A : 
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2B : 





X 
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2C : 
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X 




X 













X 













X 




X 





X 




X 





X 




X 
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3A : 
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3B : 
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3C : 
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X 




X 













X 













X 



4 



Class 4 





X 


X 


X 




X 


X 


X 




X 


X 


X 


AA : 





X 





AB : 
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AC : 


X 
















X 




X 













X 






For each class there are six possible permutations of the rows, so there are 72 individual 
cases in all [9]. 

The rate for neutrinoless double beta decay depends on the effective Majorana mass, 
which is equal to the magnitude of the v e — v e element of the neutrino mass matrix, 

|M ee | = \m lC 2 12 c 2 13 + m 2 e-^ 2 s 2 12 c 2 n + m 3 e- i<t,:i s 2 13 e 2iS \. (A) 

Note that \M ee \ can be written in a form that is independent of the Dirac CP phase S by 
redefining 3 . We choose the form above because the Majorana phases will be constrained. 
The latest experimental result from EXO-200 [TT] shows that the effective mass is less than 
140 — 380 meV at 90% C.L. In the foreseeable future, experiments such as KamLAND-Zen 
will reach a sensitivity of about 50 meV or below [12]. 

In principle, leptogenesis may provide hints of the structure of the Yukawa matrix; see 
Ref . [13] for a recent review. In what follows, we adopt a minimal scenario of leptogenesis [14] . 
We work in the single flavor approximation in which Mi > 10 12 GeV, so that the flavor 
composition of the leptons does not affect the baryon asymmetry of the universe. We also 
assume the right-handed neutrinos to be hierarchical, M 2 ,M 3 3> Mi > 10 12 GeV. Then in 
the standard model the baryon asymmetry is given by [H] 

r]m = r ^- = -9.72 x 10~ 3 x e x x ri, (5) 

n 7 

where ei is the CP asymmetry in the decay of the lightest right-handed neutrino, and rj is 
the wash out efficiency factor, which can be obtained by solving the Boltzmann equation. 
We use a very simple analytic fit [H] 

/0.55 x KT 3 eV\ 116 

■K — ^ — J • (6) 

where rh\ is 

v 2 

a=e,/i,T 1 a=e,/i,T 
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Equation © is valid for M 1 < 10 14 GeV and ih\ > 0.01 eV. With these constraints, the 
baryon asymmetry is 

, so = 1^ _ 3 .4 x 10- x (8) 

n 7 mi 



where e\ can be written as 



f > - -i^plk Im[(AA, ^ 1(AA,f)l11 



Since A = -Mi' 2 Y, we have 



Cl = -J- 9/ M * Vim (yyt)2 

3Mi \- Im(yi e ^ e + y^y*^ + ?/i r ?/* T ) 2 

I,., 12 i |„, 12 i |, # , 12 ^ U ' ) 



16^2^ |yie| 2 + IZ/l^l 2 + I?/1t 

Hence, the baryon asymmetry is proportional to the lightest right handed neutrino mass 
Mi. The observed baryon asymmetry is i]bo = 6.19 x 10~ 10 [15J. Since Y is determined by 
the light neutrino mass matrix, we can calculate the lightest right-handed neutrino mass Mi 
from Eqs. Q, ©, and (HDD, ensuring that 10 12 GeV < M x < 10 14 GeV. 



3 Phenomenology 



3.1 Class 1 (M 23 = 0) 



The mass matrices of the three textures in Class 1 always have a zero in the (2, 3) entry. In 
fact, M23 = is the only condition on Y for all the textures in Class 1, and the condition 
is the same for both mass orderings (normal and inverted). Take Class 1A for example; we 
can write 



Y 



a b 
c d 
e 



(11) 
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where a, 6, c, d, e are all nonzero complex numbers. Then the mass matrix of the light neu- 
trinos becomes 



M — Y Y 



a 2 + c 2 + e 2 ab cd 
ab b 2 

cd ci 2 



Comparing Eq. f[T2"j) with the standard parametrization, if M 2 3 
for Y as follows: 



(12) 



0, we can find a solution 



b = ±M 22 /2 , a = M 12 /b, d = ±M^ Z , c = M l3 /d, e = ±{M U - a 2 - c 



r l/2 



2 _ 2x1/2 



(13) 



Since a, 6, c, d, e are all nonzero complex numbers, a solution always exists. Solutions for Y 
for the other two textures in Class 1 may be derived in a similar fashion. Therefore, Class 1 
is defined by the necessary and sufficient condition M 2 3 = 0, which can be written as 

m 3 e l ^U T3 U^ 3 + m 2 e l(t>2 U T2 U^ 



VTL\ 



Taking the absolute square gives 



U Tl U, 



(14) 



til 



m\\U^\ 2 \U Tl \ 2 -ml\U^ 2 \ 2 \U T2 \ 2 -ml\U^\ 2 \U T3 \ 2 = 2Re(m 3 e- i ^U; 3 U; 3 m 2 e i ^U^ 2 U r2 ), (15) 
or, defining (ft = (f> 3 — (fi 2 , 

m 2 \U^\ 2 \U rl \ 2 - m 2 2 \U^ 2 \ 2 \U r2 \ 2 - m^lM^a! 2 

= -2m 2 m 3 c 2 13 s 23 c 23 Re [e~ H '(c 12 c 23 - s 12 s 23 s 13 e l5 )(c 12 s 23 + si 2 c 23 si 3 e J<5 )] 
= -2m 2 m 3 c 2 13 s 23 c 23 x 

Re [c\ 2 S 23 C 23 COS (j) + Si2Ci 2 Si 3 (c23 - S 23 ) cos(0 - 5) - s 2 2 S2 3 c 23 s 2 3 cos(0 - 25)] .(16) 



Expanding the cosines yields the form 



C = A cos (f) + B sin 0, 



(17) 



with A, B and C as listed in Table |2J 

For Eq. fTTTjl . if C 2 > A 2 + B 2 , there is no solution; if C 2 < A 2 + B 2 , there are two 
solutions: 

B ± VA 2 + B 2 -C 2 , 
2arctan ; . (18) 



A + C 
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C = A cos + B sin 



Class 


A 


B 


C 




-2m 2 m 3 Ci 3 S23C23X 


-2m 2 m3Ci 3 S23C23Sl2Sl3 x 


m?|C/ Ml | 2 |C/ Tl | 2 


1 


[C? 2 S23C 2 3 + S12C12S13(C2 3 - S 2 23 ) COS 5 


[ci 2 (4 3 - sl 3 )sin5 






-S?2 S 23C23S? 3 COs(2<5)] 


-S12S23C23S13 sin(2(5)] 


-m 2 |^ 3 | 2 |C/ r3 | 2 


2 


— 2m 2 m3cf 3 Si2C23Si 3 x 

L^i o Qno prm /[ — 1 — e, c 1 o prm (yA i 
L°12 ^ 23 ^ UB u T «->12 l '23 a 13 l-l'o^z.u 


— 2m 2 m 3 cf 3 Si2C2 3 Si 3 x 

L°12 ^ 23 u t ■JI2'^23' 5 13 


m 2 |f/ el | 2 |f/ Tl | 2 

-m 2 |f/ e2 | 2 |^2| 2 

-m 2 |f/ e3 | 2 |^3| 2 


3 


2m 2 m3C 13 Si2S23Sl3 x 

[ci 2 c 23 COS 5 - S 12 S 2 3S 13 cos(25)] 


2m 2 m 3 c 13 Si2S23 s i3 x 

[ci 2 c 23 sin 5 - S12S23S13 sin(25)] 


m 2 |f/ el | 2 |f/ Ml | 2 

-m 2 |f/ e2 | 2 |^2| 2 

-m 2 |f/ e3 | 2 |M 2 




-2m 2 1 m 3 1 cl 3 s 2 3C23 x 


-2m 2 1 m 3 1 cj 3 S23C23Sl2Sl3 x 


mr 2 |C/ M i| 2 |C/.i| 2 


4A 


[C? 2 S23C 2 3 + S 12 C 12 S 13 (4 3 - s| 3 ) COS (5 


[ci 2 {c 2 23 - s 2 23 )sin5 


-m 2 - 2 |^ 2 | 2 |^ 2 | 2 




-S?2 S 23C 2 3-S?3 C0S(25)] 


-Si 2 s 23 c 23 s 13 sin(25)] 


-m3 2 |^ 3 | 2 |f/r 3 | 2 


4B 


—zm 2 m 3 c 13 Si2C2 3 Si 3 x 

[Ci 2 S 2 3 COS 5 + Si2C 23 Sl 3 COs(2(5)] 


O/no /m ^" n^ n \/ 

—zm 2 m 3 c 13 Si2C2 3 Si 3 x 
[ci 2 s 2 3 sin 5 + si 2 c 23 si 3 sin(25)] 


mr 2 |f/ el | 2 |f/ rl | 2 

-m 2 2 \U e2 \ 2 \U r2 \ 2 

-m 3 2 \U e3 \ 2 \U T3 \ 2 


4C 


2m 2 1 m 3 1 cl 3 s 12 c 23 s 13 x 

[ci 2 c 23 cos 5 - S12S23S13 cos(2c5)] 


2m 2 1 m 3 1 cj 3 s 12 C23Si3 x 
[C12C23 sin 5 - S12S23S13 sin(2c5)] 


mr 2 |C/ e i| 2 |^i| 2 

-m 2 - 2 |f/ e2 | 2 |^2| 2 

-m 3 2 |f/ e3 | 2 |^3| 2 



Table 2: The coefficients A, B and C for each class. 



We can write Eq. (THj) in terms of a s 



mi = e 



(19) 



Since mi is a non-negative real number, we get 



'2 = ~ 



arg[ 



(20) 



and 



03 = 02 + 0- 



(21) 



Therefore, if m\ (m^) and 5 in the normal (inverted) ordering are known, we can calculate 
A, B and C using the five measured oscillation parameters in Table [TJ We scan the 5 and 
m i (^3) space to find the allowed regions, which are defined by the condition C 2 < A 2 + B 2 
(see Fig. [T]for the normal ordering and the upper panel of Fig. [2] for the inverted ordering). 
We show allowed regions corresponding to the best-fit parameters, and those allowed at 2a. 
We also plot iso-0 2 and iso-|M ee | contours using the best-fit oscillation parameters. We only 
show the contours for the plus sign of in Eq. f fT8|) because changing 5 to 360° — 5 yields 
the same contours for the minus solution. 

The allowed regions can be further constrained using leptogenesis. We assume the lightest 
right-handed neutrino has mass between 10 12 GeV and 10 13 GeV and is much lighter than 
the others so that we can use Eq. ([8]). We also require rh\ > 0.01 eV. From Eqs. ([7]) and 
( ITOj) . we see that the baryon asymmetry depends on the sign choices of in Eq. ( !T8|) but 
not on the sign choices in Eq. 013 j) . because different choices of signs in Eq. f)13p change the 
signs of all parameters in one row of the Y matrix, which yield the same baryon asymmetry. 
The baryon asymmetry also depends on the row of Y that is associated with the lightest 
right-handed neutrino mass, but the order of the other two rows does not affect the baryon 
asymmetry. 

Since Classes 1A, IB and 1C have different textures, the allowed regions for successful 
leptogenesis are also different in these three cases. Here we only consider Class 1A as an 
example. We find that successful leptogenesis is not possible for the normal ordering. For 
the inverted ordering, the allowed regions are shown in the three lower panels of Fig. |2j 
Although the constraints on 5 vary according to which right-handed neutrino is lightest, in 
no case is the lightest left-handed neutrino allowed to be above 100 meV. 
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3.2 Class 2 (M 13 = 0) 

Similar to Class 1, the only condition for Class 2 is M 13 = 0, which is independent of ordering 
and can be written as 

mi = uZu • (22) 

After taking the absolute square, then as in Class 1 this may be put in the form of Eq. (fT7|) , 
with A, B and C as in Table [2j 

The solutions for 2 and 03 then proceed as in Class 1. The allowed regions for the 
inverted ordering are shown in Fig. EJ along with iso-0 2 and iso-|M ee | contours. We see that 
the solution found in Ref. [8] with Ai e = A 2t = or Ai T = A 2e = for the inverted ordering is 
a special case of our model with 777,3 = 0. Leptogenesis predictions for Class 2A 10 are also 
shown in Fig. [31 and give an upper bound on m 3 of about 200 meV. The allowed regions for 
the normal ordering are similar to Fig. |H 

3.3 Class 3 (M u = 0) 

The only condition for Class 3 is M 12 = 0, which can be written 

mi = UJZ ' (23) 

This condition is the same for both mass orderings and as before this may be put in the 
form of Eq. (Q2D with A, B and C as in Table 

Note that Class 3 is the same as Class 2 for any ordering with U T j — > U^, or s 2 3 — > — c 2 3 
and c 23 — > S23, which is the same as 5 — > 5 + 180° when 6* 2 3 = 45°. Since # 23 ~ 45°, the 
phenomenology of Classes 2 and 3 will be similar. The allowed regions for the normal ordering 
are shown in Fig. HJ also shown are predictions for \M ee \ and </> 2 and regions compatible with 
leptogenesis. The allowed regions for the inverted ordering are similar to those for Class 2 
in Fig. [3] with 5 — > 5 + 180°. We note that the solution found in Ref. [8] with Ai e = A 2At = 
or Ai M = A 2e = for the inverted ordering is a special case in our model with m 3 = 0. 
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3.4 Class 4A ((M -1 ^ = 0) 



Y 



(24) 



Class 4 has no texture zeros in the mass matrix. However, the existence of a solution for Y 
still depends on only one condition. Take Class 4A for example, in which case 

a b c 
d 
e 

where a, b, c, d, e are all nonzero complex numbers. Then the mass matrix becomes 

a 2 ab ac 
ab b 2 + d 2 be 



M = Y Y 



ac be c 2 + e 2 



(25) 



and we see that if ^t 11 

M12 



is satisfied, the other variables can be directly determined. 



Hence the Class 4A condition is equivalent to MnM 2 3 = M12M13, which means the (2, 3) 
cofactor of M, C 23 = Mi 2 M 3 i — MnM 32 , vanishes. Since (M -1 ),;, = rrrr C 1 ^, the condition 
for Class 4 A is equivalent to (M _1 



12:', 



0, i.e., a texture zero in the inverse mass matrix. 



Since M 1 = Vdi&g(m 1 , m 2 , m 3 1 )V T ', we can write the condition as 



(26) 



or 



The allowed regions for the normal ordering are shown in Fig. |5j The allowed regions for 
the inverted ordering are similar to Fig. [1] (Class 1 NO), and the iso-</>2 and iso-|M ee | contours 
are similar with <j) 2 — > —<p2- The similarity of an 10 scenario with an NO one may seem 
unusual, but can be understood by looking at the form of the A, B, and C coefficients in 
Table 2; multiplying the coefficients for Class 4A 10 by m 2 m3, and dividing the coefficients 
for Class 1 NO by m 2 m 3 , we see that A and B are the same for the two cases. For the 
C coefficient, the dominant term in each case is the third one, proportional to |t^3| 2 |^ T 3| 2 
times the ratio of a larger mass to a smaller one. 

The comparison of Class 4A NO with Class 1 10 is more nuanced: when the lightest 
mass (mi for NO, m 3 for 10) is very small, the first two terms in the C coefficient have 



m- 
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similar size for Class 1 10, but only the first term is dominant for Class 4A NO. However, 
when the lightest mass is not so small, such that m\ ~ m 2 in the NO, then the same terms 
in the C coefficient are dominant. Thus for higher values of the lightest mass, although not 
necessarily so large that all three masses are quasi-degenerate, Classes 4A NO and 1 10 give 
similar predictions. This can be seen by comparing Figs. 5 and 2: although the allowed 
regions and contours are quite different when the lightest mass is below 20 meV, note the 
similarity of the \M ee \ = 100 meV and 2 = 60° contours. 

3.5 Class 4B ((M _1 )i 3 = 0) 

Similar to Class 4A, the condition for Class 4B is (M _1 ) 13 = 0, which can be written as 



u Tl u el 

This condition is the same for both mass orderings, and the analysis follows as in previous 
classes, with A, B and C given in Table |2j 

The allowed regions for the normal (inverted) ordering are shown in Fig. [6] (Fig. [7]). The 
inverted ordering for this case is also similar to Class 2 NO: multiplying the A, B, and C 
coefficients by m 2 7«3 for Class 4B IO and dividing them by m 2 m 3 for Class 2 NO, A and B 
are identical for the two cases, and the dominant terms in C are also the same. As was true 
in the previous section, the reverse correspondence between Class 4B NO and Class 2 IO 
exists only for larger values of the lightest mass (see Figs. 6 and 3 when the lightest mass is 
above 50 meV). 

3.6 Class 4C ((M _1 )i 2 = 0) 

Similar to Class 4A, the condition for Class 4C is (M" 1 )i2 = 0, which can be written as 

The corresponding values of A, B and C in Eq. (TTTjl are given in Table |2j 

Note that Class 4C is the same as Class 4B with U T j — > U^, or s 23 —> — c 23 and c 2 3 — > s 2 3- 
As noted in Sec. 3.3, this transformation is equivalent to 5 — > 5 + 180° when 823 = 45°. 
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Therefore the allowed regions of Class 4C are similar to Class 4B in Fig. [6] for the normal 
ordering and in Fig. [7] for the inverted ordering with 5 — > 5 + 180°. 

The inverted ordering for this case is also similar to Class 3 NO, as can be seen by 
examining the A, B and C coefficients, and Class 4C NO and Class 3 10 give similar results 
for larger values of the lightest mass. Thus there is a general pattern that the texture zero 
NO and corresponding cofactor zero 10 have similar predictions, and texture zero 10 and 
corresponding cofactor zero NO have similar predictions when the lightest mass is not too 
small. 

4 Conclusions 

We extended the most economical type I seesaw model to include three right-handed neu- 
trinos. The simplest cases that fit the data have four texture zeros in the Yukawa couplings 
that connect the left-handed and right-handed neutrinos. These are equivalent to a single 
texture or cofactor zero for an off-diagonal element of the light neutrino mass matrix M. 
The cofactor zero condition is itself equivalent to a texture zero in M _1 . We used the latest 
experimental data to obtain the allowed regions for the lightest neutrino mass and Dirac CP 
phase 5, which can be measured in future neutrino experiments. We also used leptogenesis 
to further constrain the allowed regions; in general there is an upper bound on the lightest 
neutrino mass of about 100-200 meV for a single-flavored leptogenesis scenario. Figures 2 
to 7 show that in any given model, not all values of 5 are consistent with the leptogenesis 
predictions. Therefore a precise measurement of 5 will reduce the number of viable models. 

Once the lightest neutrino mass and Dirac CP phase are determined, these models make 
definite predictions for neutrinoless double beta decay. From the iso-|M ee | contours in Figs. 
1-7 we see that \M ee \ is generally proportional to the lightest mass. This behavior is clearly 
evident for the quasi-degenerate spectrum. However, for some classes, \M ee \ is strongly 
dependent on the Dirac CP phase 5 and is not determined by the lightest mass alone. We 
plot the allowed regions in the (|M ee |, lightest mass) plane for these classes separately; see 
Figs. El EI US and CU 

The minimum value of \M ee \ for the best-fit oscillation parameters and the 2a lower 
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Table 3: The minimum values of \M ee \ (in 10 3 eV) in each class for the best-fit oscillation 
parameters, and the 2a lower bounds. 



Class 


Best-fit 


2cr lower bound 


NO 


IO 


NO 


IO 


1 


142.5 


19.0 


129.8 


15.4 


2 


1.4 


46.9 


0.3 


44.8 


3 


0.0 


47.4 


0.0 


45.2 


4A 


0.0 


150.3 


0.0 


138.3 


4B 


6.1 


18.4 


5.7 


15.1 


4C 


8.4 


18.2 


7.5 


14.9 



bounds are shown in Table 3. We find that the Class 1 NO and Class 4 A IO have a 
minimum value for \M ee \ of around 150 (130) meV for the best-fit parameters (at 2a), and 
could therefore be excluded by the 0vj3j3 decay experiments that are currently running; 
Classes 2 IO and 3 IO have a minimum \M ee \ of around 50 meV and can be definitively 
tested by experiments under construction. For Classes 3 NO and 4A NO the current lower 
bound on \M ee \ is zero, and given the current measurements of the oscillation parameters, 
0vf3f3 decay will not constrain them. The remaining models have a minimum \M ee \ in the 
range 1 — 20 meV and can only be completely probed by significant improvements in the 
sensitivity of 0u/3/3 experiments. The sum of light neutrino masses can also be used to 
provide additional discrimination among these models. However, there is a general pattern 
that a texture zero NO and the corresponding cofactor zero IO have similar predictions, and 
a texture zero IO and the corresponding cofactor zero NO have similar predictions when the 
lightest mass is not too small. Therefore it may be difficult to uniquely specify the model 
from data, although experiments designed to determine the mass ordering could resolve this 
ambiguity. 

Since the models studied in this paper are equivalent to a single texture or cofactor 
zero for an off-diagonal element of M„, one might also consider examining models with a 
single texture or cofactor zero for a diagonal element of M v . Although not obtainable from 

14 



texture zeros in the Yukawa couplings in a type I seesaw model, such models also have seven 
parameters in the light neutrino mass matrix and can be analyzed in a similar fashion. We 
will study the phenomenology of and possible motivation for these models in a subsequent 
paper. 
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Figure 1: The allowed regions in the (mi, S) plane for Class 1 and the normal ordering. 
The dark shaded regions correspond to the best-fit parameters of the oscillation parameters, 
while the light shaded regions are allowed at 2a. The solid lines are iso-|M ee | contours (in 
meV) and the dashed lines are iso-0 2 contours (in degrees). 
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m 3 (eV) 

Figure 2: The upper panel shows the allowed regions in the (m.3, 5) plane for Class 1 and the 
inverted ordering. The shading and line types in the upper panel are as in Fig. 1. The lower 
panels show the allowed regions for Class 1A and the inverted ordering with the additional 
constraint of successful single-flavored leptogenesis. The hatched (dark shaded) regions use 
the plus (minus) solution of <p2- From left to right the three graphs have the first, second 
and third row of Y associated with the lightest right-handed neutrino mass, respectively. 
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Figure 3: Same as Fig. [21 except for Class 2 and 2A and the inverted ordering. 
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Figure 4: Same as Fig. [2j except for Class 3 and 3A and the normal ordering. 
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Figure 5: Same as Fig. |2l except for Class 4A and the normal ordering. 
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Figure 6: Same as Fig. |2j except for Class 4B and the normal ordering. 
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Figure 7: Same as Fig. [21 except for Class 4B and the inverted ordering. 
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Figure 8: The allowed regions in the (\M ee \,m 3 ) plane for Class 2 and the inverted ordering. 
The dark shaded regions correspond to the best-fit parameters of the oscillation parameters, 
while the light shaded regions are allowed at 2a. 
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Figure 9: Same as Fig. [HI except for Class 3 and the normal ordering. 
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Figure 10: Same as Fig. [HI except for Class 4A and the normal ordering. 
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Figure 11: Same as Fig. EJ except for Class 4B and the inverted ordering. 
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